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Abstract. This report is an update and extension of our paper accepted for publication in ApJ ( arXiv:0802.4262 ). Since 
objects at the same redshift should have the same luminosity distance and the distance moduli of type la supernovae (SNe 
la) obtained directly from observations are completely cosmology independent, we obtain the distance modulus of a gamma- 
ray burst (GRB) at a given redshift by interpolating or iterating from the Hubble diagram of SNe la. Then we calibrate five 
GRB relations without assuming a particular cosmological model, from different regression methods, and construct the GRB 
Hubble diagram to constrain cosmological parameters. Based upon these relations we list the cosmology-independent distance 
moduli of 42 GRBs between redshift of 1.44 and 6.60, with the 1-C7 uncertainties of 1-3%. 
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INTRODUCTION 

Gamma-ray burst (GRB) luminosity/energy relations are connections between measurable properties of the prompt 
gamma-ray emission with the luminosity or energy. In recent years, several empirical GRB luminosity relations as 
standard candles for cosmology research at very high redshift have been proposed Hlfl, |2(]. An important concern in 
the application of GRBs to cosmology is the dependence on the cosmological model in the calibration of GRB relations 
in many works. For the difficulty to calibrate the relations with a low-redshift sample, GRB relations have usually been 
calibrated by assuming a particular cosmological model. Therefore the circularity problem can not be avoided easily. 
Many previous works treated the circularity problem by means of statistical approaches. However, we note that the 
circularity problem is not circumvented completely by means of statistical approaches, because a particular cosmology 
model is required in doing the joint fitting. 

In this report, we present two new methods (interpolation method and iterative method) to calibrate the GRB 
relations in a cosmological model-independent way. This report is an update and extension of our paper accepted 
for publication in ApJ|3]. It is obvious that objects at the same redshift should have the same luminosity distance in 
any cosmology. There are so many SNe la that we can obtain the luminosity distance (also the distance moduli) at 
any redshift in the redshift range of SNe la by interpolating from SN la data. Recently, on the basis of smoothing the 
noise of supernova data over redshift, the authors in (01, iS] suggested a non-parametric method in a model independent 
manner to reconstruct the luminosity distance at any redshift in the redshift range of SNe la by the iterative method. 
Furthermore, the luminosity distance of SNe la obtained directly from observations are completely cosmological 
model independent. Therefore, we can obtain the distance moduli of GRBs in the redshift range of SNe la and calibrate 
GRB relations in a completely cosmological model independent way and use the standard Hubble diagram method to 
constrain the cosmological parameters from GRB data at high redshift obtained by utilizing the relations. 



COSMOLOGICAL MODEL INDEPENDENT HIGH-Z GRB DISTANCE MODULI 

We first calibrate five GRB luminosity/energy relations with the sample at z < 1.4, i.e., the luminosity (L)-spectral 
lag (Tiag) relation |6], the L-variability (V) relation [7], the L-E^ relation |8], the collimation-corrected energy (Ey)-Ep 
relation [9], the Trt-^ relation |2], where Trt is the minimum rise time in the GRB light curve. A GRB luminosity 
relation can be generally written in the form of 



logy = a + blogx, 



(1) 



FIGURE 1. The Hubble Diagram of 192 SNe la 
(red dots) and the 69 GRBs obtained by using the 
cosmology-independent methods. The 27 GRBs at 
z < 1.4 are obtained by using the interpolation and 
iterative methods directly from SN la data; and the 
42 GRBs at z > 1.4 are obtained by utilizing the five 
relations calibrated with the sample at z < 1.4 us- 
ing by the cosmology-independent methods (black 
circles: the interpolation method; blue stars: the it- 
erative method). The curve is the theoretical dis- 
tance modulus in the concordance model (wo = — 1, 
O.M = 0.27, fl^ = 0.73), and the vertical dotted line 
represents z= 1 .4. 

where a and b are the intercept and slope of the relation respectively; y is the luminosity (L in units of erg s^') 
or the collimation-corrected energy (Ey in units of erg); x is the GRB parameters measured in the rest frame, e.g., 
Tiag(l +z)"V(0-l s), y(l +z)/0.02,£p(l +z)/(300keV), Trt(1 +z)-V(0.1 s). (We adopt the data of these quantities 
fromRef. |2].) 

We adopt the data of 192 SNe la ifioll and show them in Figure 1. There is only one SN la point (the redshift of 
SN1997ff is z = 1.755) at z > 1.4, therefore we exclude it from our SN la sample used to interpolate the distance 
moduli of GRBs in the redshift range of SN la sample. 

Following a well known procedure in the analysis of large scale structure, Shafieloo et al |Q] used a Gaussian 
smoothing function rather than the top hat smoothing function to smooth the noise of the Sne la data directly. In 
order to obtain important information on interesting cosmological parameters expediently, when doing the Gaussian 
smoothing|4] lnt/i(z), rather than the luminosity distance di{z) or distance modulus ji{z), is studied by the iterative 
method. We thus follow the iterative procedure and adopt the results from Ref. |5], 

lndL{zY„=lndLizyi,_, +A^(z)i:(lnjf^-fe) -ln«'L(z,)Li)exp [- (W ((1 +z)/(l +z,)))/(2A2)] , (2) 

i 

with a normalization parameter A^(z)^^ = L/Cxp [ — (in^ ((1 + z)/(l + /(2A^)] , and A = 0.6. diizYl represents 
the smoothed luminosity distance at any redshift z after the nth iteration and dL{z)Q denotes a guess background 
model. It has been shown that the results are not sensitive to the chosen vaule of A and the assumed initial guess 
model. df"^{zi) is the observed one from the SN la data. The best fitting result is obtained by minimizing = 

The isotropic luminosity of a burst is calculated by L = A%dll\,o\o, where c/l is the luminosity distance of the burst 
and fboio is the bolometric flux of gamma-rays in the burst. The isotropic energy released from a burst is given by 
£iso = 47rc/L'5boio(l where 5boio is the bolometric fluence of gamma-rays in the burst at redshift z- The total 

coUimation-corrected energy is then calculated by Ey = Fbeam^^iso, where the beaming factor, Fbeam is (1 — cos 0jet) 
with the jet opening angle (0jet), which is related to the break time (fb). 

We determine the values of the intercept (a) and the slope ib) calibrated with the GRB sample at z < 1 .4 by using the 
interpolation! 3 1 and iterative Q methods, respectively. We first use the same regression method as used in Ref.|01: the 
bisector of the two ordinary least-squares(OLS) linear regressions, which has been discussed in Ref. [|1 111 , OLS(YIX) 
(OLS regression of the dependent variable Y against the independent variable X) and its inverse OLS(XIY). In order to 
avoid specifying "dependent" and "independent" variables, the OLS(YIX) and the OLS(XI Y) lines should be bisected. 
The OLS regressions does not take the errors into account; but the use of weighted least-squares(WLS), taking into 
account the measurable uncertainties, results in almost identical best fits. Therefore, the WLS bisector regression can 
be obtained by bisecting the WLS(YIX) and the WLS(XIY) lines. The calibration results from two regression methods 
(the OLS bisector and the WLS bisector) are summarized in Table 1 and we also list the results calibrated with the 
same sample (z < 1 .4) by assuming the ACDM model for comparison. From Table 1, we find that results from the two 
regression methods make no significant difference from each other Therefore, taking into account the measurement 




TABLE 1. Calibration results (a: intercept, b: slope) with their 1-C7 uncertainties for the five GRB luminos- 
ity/energy relations with the sample at z < 1.4, by using the interpolation and iterative methods directly from 
SNe la data, and by assuming a particular cosmological models (the ACDM model). Only the data available 
from Table 4 in Ref.|2] are used to calibrate the five GRB relations. For each relation, the results on the first 
and second lines are calibrated from the OLS and WLS bisector regressions, respectively. 



Interpolation method 



Iterative method 



ACDM model 



Xiag-L relation 


52.22±0.09 -1.07±0.13 
52.22±0.09 -1.07±0.14 


52.22±0.09 -1.10±0.14 
52.22±0.09 -1.10±0.15 


52.15±0.10 -1.11±0.14 
52.14±0.10 -1.10±0.14 


V-L relation 


52.59±0.13 2.05±0.27 
52.59±0.13 2.08±0.37 


52.56±0.13 2.03±0.30 
52.57±0.13 2.10±0.42 


52.48±0.13 2.04±0.30 
52.49±0.13 2.11±0.42 


L-Ep relation 


52.26±0.09 1.69±0.11 
52.26±0.09 1.69±0.16 


52.23±0.09 1.65±0.12 
52.23±0.09 1.65±0.16 


52.15±0.09 1.66±0.12 
52.15±0.09 1.65±0.16 


Ey-Ep relation 


50.71±0.07 1.79±0.18 
50.71±0.06 1.68±0.29 


50.68±0.06 1.75±0.18 
50.68±0.06 1.63±0.30 


50.59±0.06 1.75±0.19 
50.60±0.06 1.64±0.30 


Trj-L relation 


52.64±0.11 -1.32±0.15 
52.64±0.11 -1.33±0.16 


52.60±0.10 -1.30±0.15 
52.60±0.10 -1.31±0.16 


52.52±0.11 -1.30±0.15 
52.52±0.11 -1.31±0.17 



uncertainties in the regression, indeed will not change the fitting parameters significantly, when the measurement 
uncertainties are smaller than the intrinsic error [2]. We also find that results obtained by assuming the ACDM with 
the same sample differ only slightly from, but still fully consistent with those calibrated by using our interpolation 
method. The reason for this is easy to understand, since the ACDM is fully compatible with SN la data. Nevertheless, 
it should be noticed that the calibration results obtained by using the interpolation or iterative method directly from 
SN la data are completely cosmology independent. 

By utilizing the calibrated relations at high redshift (z > 1 .4)^e are able to obtain the luminosity (L) or energy (Ey) 
of each burst at z > 1 .4. We use the same method used in Ref.[2] to obtain the best estimate ji for each GRB which is 
the weighted average of all available distance moduli. The derived distance modulus for each GRB is 

M = (LM./cT^,)/(i:(T^f), (3) 

i i 

with its uncertainty — <7^^) ''^^^ where the summations run from 1 to 5 over the five relations used in Schaefer 
(2007) with available data. 

We plot the Hubble diagram of the 69 GRBs obtained by using the interpolation and iterative methods in Figure 
1 . The 27 GRBs at z < 1 .4 are obtained by using the interpolation and iterative methods directly from SNe data. The 
distance moduli of these 42 GRBs at z > 1 .4 are obtained by utilizing the five relations calibrated with the sample 
at z < 1.4 using the cosmology-independent methods. The derived distance moduli for the 42 GRBs (z > 1.4) with 
different methods are listed in Table 2, together with the average values between different methods, which should be 
the least biased and most robust values to be used to study the expansion history of the Universe up to z = 6.60. It 
should be noted that the \-o uncertainties listed in Table 2, ranging between 1% to 3%, include both the measurement 
uncertainties and intrinsic scattering in these luminosity/energy relations. 

In Figure 2 we show examples of cosmological parameter fitting by the minimum method. Figure 2a and b show 
the joint confidence regions for (ilMi^^A) in the ACDM model from the distance moduli of these 42 GRBs (z > 1.4) 
obtained by utilizing the five relations calibrated with the sample at z < 1.4, using the interpolation and iterative 
methods, respectively. Figure 2c and d represent the contours of likelihood in the (i2M,wo) plane in the dark energy 
model with a constant wq for a flat universe. Here we adopt Hq = 70 km s^^Mpc^^. All fitted parameters are listed 
in Table 3. We find that the fitting results from the OLS regression differ only slightly from, but still fully consistent 
with those from the WLS regression. But the use of WLS, taking into account the measurement uncertainties in the 
regression, results in almost identical best fits. 



TABLE 2. Cosmology-independent distance moduli of 42 GRB (1.44 < z < 6.60). jUi and 1x2 are 
obtained by the interpolation and iterative methods, respectively; fx" and fi^ are obtained from the OLS 
and WLS bisector regressions, respectively. ^ = [n^ + fxf + + fi^) /4 and (J^ = (ct^« + a^t + 0^ + 
) /4, which should be least biased and most robust for applications. The 1-CT uncertainties include both 
the measurement uncertainties and intrinsic scattering in these luminosity/energy relations. 



GRB z /xf M 
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.75 


44, 


.62 


± 
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.57 


± 


0, 
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SUMMARY 

Since the distance modulus of any SN la is completely cosmological model independent, we can obtain the distance 
modulus of a GRB at a given redshift by interpolating or iterating from the Hubble diagram of SNe la at z < 1.4, in 
order to calibrate the GRB luminosity relations in a completely cosmology independent way. Since our method does 
not depend on a particular cosmological model when we calibrate the parameters of GRB luminosity relations, the 
so-called circularity problem can be completely avoided. With this method, we obtained the cosmology-independent 
distance moduli of 42 GRBs between redshift of 1 .44 and 6.60, which are listed in Table 2. 
With these cosmology-independent distance moduli of high redshift GRBs, we construct the GRB Hubble diagram 
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FIGURE 2. (a) and (b) Joint confidence regions for (flM,fiA) in the ACDM model from tiie data for 42 GRBs (z > 1.4) obtained 
by utilizing the five relations calibrated with the sample at z < 1 .4, using the interpolation and iterative methods, respectively. The 
dashed line represents the flat universe, (c) and {d) Contours of likelihood in the (f2M:>i'o) plane in the dark energy model with 
a constant wq for a flat universe. For each panel, the horizontal line is for a prior of £ly^ = 0.27. For all four panels, the plus 
sign indicates the best-fit values and the contours correspond to 1,2, and 3-(T confidence regions. The numerical values of fitted 
parameters are listed in Table 3. 



TABLE 3. The cosmological fitting results from the distance moduli of these 42 GRBs (z > 1.4) 
obtained utilizing the five relations calibrated with the sample at z < 1.4 using the interpolation and 
iterative methods, from different linear regressions. The fitted is for the flat ACDM model. For the 
dark energy model with a constant equation of state vvq is for a prior of £ly^ = 0.27. For the every two 
lines of each parameter, the results on the first line are derived from the interpolation method, and that 
on the second line are from the iterative method. 





OLS(YIX) 


WLS(YIX) 


OLS(XIY) 


WLS(XIY) 


OLS Bisector 


WLS Bisector 




35+" '"' 

"■-'■-'-0.05 

n ^s+o ('7 

0.38_„()7 


35+" '^'' 
"•■-'■'-0.06 

n ^s+0.07 

U.38_oo5 


0.1S_,)„7 
"■^^-0.08 


21 +"■"'' 
"•^^-0.07 

"•^^-0.08 


"•^'-0.06 

30+''^'"' 

"•■'"-0.06 


n 7S+0 05 
"■2*^-0.06 

31+°'"' 

"■■'^-0.06 


WO 


-0 72+"^2i 
"•'^-0.30 

-0 63+°^2i 


-0 73+'''23 
"•'-'-0.31 

-0 63+'^-2'' 
"■''-'-0.26 


-1 51 +'^■5" 
^•^^-1.20 

-1 25+''^3'' 

^•'^-'-0.75 


1 00+0.43 
^•^^-0.96 

-1 09+''-^5 
^■"^-0.66 


-0.97^°i« 
-0.84^°;^4 


-0.94l°:^v 
-0.82^0i^ 



and constrain cosmological parameters by the minimum method as in SN la cosmology. We obtain Hm = 0.28+Q Qg, 
Ha = 0.72+Q Q5 for the flat ACDM model from the GRB data obtained by using the interpolation method, and 
= 0.31+Q Qg, Ha = 0.69+Q Qg from the data obtained by using the iterative method. For the dark energy model 
with a constant equation of state, ■we obtain wq = — 0.94+Q4g and wo = — 0.82+Q3g for a flat universe from the data 
obtained by the two methods respectively, which is consistent with the concordance model within the statistical error. 
Our result suggests the the concordance model (wq = — 1, iijvi = 0.27, Ha = 0.73) is still consistent with the GRB 
data at higher redshift up to z 6.6. 
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